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Abstract 



^ We carry out further analysis of the Hamiltonian approach to Yang-Mills theory in 2+1 

dimensions which helps to place the calculation of the vacuum wave function and the string 
tension in the context of a systematic expansion scheme. The solution of the Schrodinger 
equation is carried out recursively. The computation of correlators is re-expressed in terms 
of a two-dimensional chiral boson theory. The effective action for this theory is calculated 
to first order in our expansion scheme and to the fourth order in a kinematic expansion 
parameter. The resulting corrections to the string tension are shown to be very small, in 
the range —0.3% to —2.8%, moving our prediction closer to the recent lattice estimates. 



1 Introduction 

A few years ago, we initiated a Hamiltonian approach to gauge theories in (2 + 1) dimen- 
sions pQ. In the Aq = gauge, which is appropriate for a Hamiltonian analysis, the complex 
components of the spatial gauge field, viz., A z , A Z: were parametrized as A z — — Z M M , 
A z = M^ l d z M\ where M, Ml are SL(N, C)-matrices for an SU{N)-gauge theory. The 
hermitian matrix H = M'M then gives the gauge-invariant degrees of freedom. The Jaco- 
bian for the change of variables from (A z , A z ) to H was explicitly calculated. This also led 
to the computation of the volume element on the physical configuration space, and hence 
the inner product of wavefunctions, in terms of the WZW action for the field H. The 
Hamiltonian was then obtained in terms of the current J = (N/ir)d z H H -1 of the WZW 
action, and the vacuum wave function \&o was calculated from the Schrodinger equation up 
to terms which are quadratic in the current J in log^o- The vacuum expectation value of 
the Wilson loop operator could then be evaluated using this wave functional. For a Wilson 
loop in the representation R, the result was 



where e is the coupling constant, and cr, ca denote the quadratic Casimir values for the 
representation R and for the adjoint representation, respectively. This value of the string 
tension is in very good agreement with lattice estimates O HI |5] . Our Hamiltonian analysis 
has also been extended to the Yang-Mills-Chern-Simons theory; the wave functions helped 
to clarify some issues regarding the dynamical contribution to the gauge boson mass and 
screening of fields [6] . 

Some of the more recent results in this approach include: 

1. A proposal for the wave function, which is very close to ours, but somewhat different, 
was made by Leigh, Minic and Yelnikov [7]. Using this wave function, an estimate 
of glueball masses was made. The values come out to be close to the lattice results 
with differences of a few percent. (It should be noted that the lattice values also have 
significant errors, especially for the higher glueballs.) 

2. The Hamiltonian formalism was extended to include scalar fields [8j. The screening 
of Wilson loops in the adjoint and other screenable representations can be related 



(W R (C)) = exp[-<rnAc] 



(1) 



where Ac is the area of the loop C. The string tension or was obtained as [2] 




(2) 
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to the formation of a bound state between a heavy scalar field and a light degree of 
freedom (the glue part). The energy of this glue-lump bound state (related to the 
string-breaking point) was calculated and shown to agree with lattice estimates to 
within ~ 9%. 

3. The formalism was developed for Yang-Mills theory on i x 5 2 [9 , motivated by 
the possibility of connection to gravity-gauge duality |10j . and also as a first step in 
developing the method for the torus, i.e., for IxS 1 X S 1 . The results for the torus 
can be useful for understanding finite temperature effects and deconfinement. 

In the light of these results, it is important to formulate a systematic expansion for the 
wave function and develop a calculational scheme for corrections to the string tension. This 
is the subject of the present paper. The recent lattice calculation of Bringoltz and Teper 
show that the string tension differs from the prediction of Q at large ./V only by about 
0.98% — 1.2% gj. Nevertheless, the deviation is statistically significant. A different lattice 
method gives a value of string tension which differs from ^ by about 1.55% [5]; again the 
deviation is considered statistically significant. These results provide another motivation 
for this paper. 

In terms of an analytical computation, there are two types of corrections to the string 
tension which may be exemplified by the two graphs shown below, where the rectangle 
of solid lines denotes the Wilson loop. The wavy lines represent the propagator or two- 
point function for the current J. The first diagram may be viewed as a correction to the 
propagator. This correction may be evaluated by computing the corrected propagator in 
the sense of an effective field theory. (The corrected propagator can then be used for the 
evaluation of the expectation value of Wilson line.) As a result, it is independent of the 




Figure 1: Examples of representation-independent (left) and representation-dependent 
(right) corrections to string tension 
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representation of the Wilson line. The second diagram is a correction which cannot be 
evaluated independently of the Wilson line. It will therefore depend on the representation 
of the Wilson line. Both types of corrections are important. The second type of corrections 
can be important in understanding screening of Wilson lines and string breaking. Some 
aspects of this problem have been treated elsewhere [8|, although a full analysis remains 
difficult. 

In this paper, we will focus on corrections of the first type, namely, those which can be 
understood as corrections to the propagator. Since the representation-dependent corrections 
are not considered, the results should be interpreted as applying in the absence of string- 
breaking, i.e., when all representations have confining area law. This means that any 
comparison with lattice value is best done at large N. 

We start with a short review of the Hamiltonian approach in section 2. Section 3 gives 
the basic expansion scheme and the procedure for calculating higher order corrections. In 
section 4, we discuss the calculation of the corrections to the string tension. A conceptual 
subtle point is discussed in appendix A and the details of the computation of the vari- 
ous corrections are given in appendix B. The paper concludes with a short summary and 
discussion. 



2 A short review 



As is standard in a Hamiltonian analysis, we use the Aq = gauge. The spatial components 
of the gauge potential can be combined as A = A z = \{A\ + 1A2), A = A z = \{A\ — 1A2). 
These can be parametrized as 

A = -dMM~ 1 , A = M t-1 5M t . (3) 

M is a complex matrix which is an element of G , the complexification of G which is the Lie 
group in which the gauge transformations take values. We will be considering G = SU(N), 
so that G c = SL(N,C). 

Time-independent gauge transformations act on M via , 

M{x) -» g{x)M(x) , g(x) e SU(N) (4) 

so that the hermitian matrix H = M^M is gauge-invariant. 

The Jacobian of the tgransformation A,A^H can be explicitly evaluated; the volume 
element for the gauge-invariant configuration space is givenby 

dfi(C)= [ d^{H) exp(2c A S wzw (H)) (5) 
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where S WZW {H) is the Wess-Zumino-Witten action for the hermitian field H and dfi(H) is 
the Haar measure for H viewed as an element of SL(N,C)/SU(N). The WZW action is 
given by 

S WZW {H) = ~J Ti(dHdH- 1 ) + ~j ^ a Ti{H- l d, L HH- l d u HH- l d a H) (6) 

Wave functions are gauge-invariant and have an inner product defined by square- 
integrability using §5§. 

The kinetic energy operator can be written in terms of its action on functionals of H as 

e 2 r S 2 e 2 



T = ~-J = ~ J u n rs (u,v)Pr(u)Ps(v) (7) 

n rs (u,v) = / g ar (x,u)K ab (x)G bs (x,v) 

J X 

where K a b = 2TY{t a Ht H~ 1 ) is the adjoint representative of H, and p a and p a are left and 
right translation operators for defined by 

[ Pa (x),M(y)} = M(y)(-it a )5^(x-y) 

[p a (x),M\y)\ = (- ita )M\y)5^(x-y) (8) 

These may be considered as functional differential operators. The Green's functions in ([7]) 
are given by 

Gma&y) = y) [S ma ~ e~\^ 2 ^(K(x,y)K~\y,y)) n 

G ma (x,y) = \ [5 ma -e-^ 2 ^{K- 1 (y,x)K(y,y)) ma ] (9) 

These are the regularized versions of the corresponding Green's functions 

G(x, y) = 1 , G(x, y) = } (10) 
vr(x - y) tt(x - y) 

The parameter controlling the regularization, e, acts as a short-distance cut-off. Expression 
([7]) is to be used on functionals where the point-separation of various factors is much larger 
than y^. 

As discussed in [T|, (M, M+) and {MV(z}, V(z)M^) give the same gauge potentials 
A, A, where V, V are, respectively, antiholomorphic and holomorphic in the complex 
coordinates z = x\ + 1x2 and z = x\ — 1x2- To avoid this ambiguity of parametrization, 
physical observables in the theory should therefore satisfy the holomorphic invariance 

H -► V{z)HV{z) (11) 
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The regularization used in Q respects this invariance. 

Wave functions can be taken to be functionals of the current J = (ca/^OHH" 1 , where 
ca is the quadratic Casimir invariant defined by CA& ah = f amn f bmn - : CA = N for SU(N). 
The action of T on wave functions of the form ty( J) can be expressed as 



e 2 c A 
2tt 



SJ a (z) 



+ 



n ab (zw)— — 



*(J) (12) 



where 



A ra (w,z) = -[d z Il rs (w,z)} Kj(z) 

with V w a b = ^-d w 5 a b + ifabcJc(w). For small e, T can be further simplified as 

5 5 



(13) 



T YM ^{J) = m 



where m = e 2 CA/27r. 

The potential energy can be written in terms of the current as 



*(J) + 0(e) (14) 



Vym 



7T 



mcA 



d 2 x : dJ a {x)dJ a {x) : 



(15) 



In terms of a regularized form with normal ordering, this expression is to be interpreted as 

c^dimG 



Vym 
<j(x, y; A) 



7T 

mcA 
1 



a(x, y; X)dJ a (x)(K(x, y)K (y, y)) a bdJ b (y) 



7T 2 X 2 



7rA 



exp 



\x- y\' 
A 



(16) 



where cr(x, y; A) is a regularized ^-function, A is the parameter of regularization, and we 
should take the limit where y/X <C 1/e 2 . The operator U ab (x,y) = [K(x,y)K~ 1 (y,y)] ab is 
such that the regularized expression for Vym satisfies holomorphic invariance. 

The total Hamiltonian, which we shall use in what follows, is thus given by 

1 5 5 



H 



in 



6 mcA 



TT 



mcA 



8J a {z) tt 2 J z w (z - w) 2 8J a {w) SJ a (z) 
d 2 x : dr(x)dJ a (x) : +im f f abc ^-^r-^ < 17 > 

Jz,w 7T{Z - W) dJa{w) dJ b (z) 
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The Schrodinger equation for the vacuum wave function was analyzed in [TJ |2] and the 

! q 

leading term in a strong-coupling expansion was obtained as = e 2 > where 



4vr 2 

S(H) = I 8J a 



m + \/m 2 — V 2 ) 
2f abc I f {3) (x,y,z)J a (x)J b (y)J c (z) + 0(J 4 ) (18) 



The function f( 3 \x,y, z) is given in [2]. 

Notice that if one restricts to modes of J with momentum <C e 2 , 

a. a to. 



2ir 2 r - 

s{H) « — / Bj a dJ 

me z c A J 



A 

Computation of expectation values reduces in this limit to a calculation in a Euclidean 
two-dimensional Yang-Mills theory with a coupling g 2 = me 2 . This was the approximation 
used in arriving at the formula ^ for the string tension. 

3 A systematic expansion scheme 

We will start with an outline of our method for calculating corrections to the formula for 
the string tensions given in ([2]). We shall first rewrite the derivation of the vacuum wave 
function as a recursive procedure for the solution of the Schrodinger equation which will 



make it clear that ( 18 ) is the lowest order result in a systematic expansion. For this purpose, 
after a rescaling of the current J, we will treat m and e as independent parameters, setting 
m = e 2 c J 4/2-7r only at the end of all calculations. In terms of these parameters, = eP > 
where P is a power series in e. (However, this is still quite different from perturbation 
theory since m is included exactly in the lowest order result for P. This recursive procedure 
is something like a resummed perturbation theory. The resummation involves collecting 
A, A in an appropriate series to define J and then including m at the lowest order which is 
another series.) 

The calculation of averages will involve integration of ^q^o over au field space. Since 
J is not an independent variable, we will express the integration in terms of a chiral boson 
field; this transformation is analogous to the fermionization of the WZW model. We can 
calculate the corrections to the J 2 -term (and eventually the string tension) by viewing this 
version in terms of the chiral boson field as a two-dimensional field theory. There will be two 
sets of contributions which are corrections to the J 2 -term to any given order in e. One term 



7 



will be a direct contribution from the recursive procedure. The other set of terms will be 
loop corrections of the two-dimensional field theory. These latter terms can be considered 
as Feynman diagrams. Vertices corresponding to currents in these loops will carry powers 
of m/Ek = m/Vk 2 + to 2 , where k is some typical momentum. These factors suppress the 
contribution of the loop integral, and so, it is advantageous to group loop corrections by 
the number of powers of m/ E^. Our procedure thus involves three steps: 

1. Solving the Schrddinger equation as a power series in e, after rescaling J, and treating 
to and e as independent parameters; 

2. Evaluating the loop corrections in the two-dimensional field theory used for computing 
expectation values; 

3. Grouping loop corrections by powers of m/Ek and calculating all the contributions 
for each power of m/Ek. 

3.1 The recursive solution 

We now turn to the first step in the expansion scheme outlined above, namely, the recursive 
solution of the Schrddinger equation. For this, we shall rescale the current as J — > 



The Hamiltonian (17) now takes the form 



/§ 2 f 1 

Ja{z) +- 
5J a (z) TT J u 



dJ a (z) TT J w z (Z - w) 2 5J a (w) 5J a (z) 

+ l - j : Bj a (z)dJ a (z) : (20) 



r,i +,v / ./;,„- JC{w) 



w,z 



F 



'tx{z — w) 5J a {w) 5Jb{z) 
The vacuum wave function is taken to be of the form = exp(^F), where 

J /il(x!,X 2 ) J^{x X )J^{x 2 ) + ^ f£l a3 ( Xl ,X 2 ,X 3 ) J^{x X )J^{x 2 )J^{x 3 ) 

+ J f£la,a<(.xi,X2,X3,xt) J^(x 1 )J a ^x 2 )J a ^x 3 )J a ^x i ) + ••• (21) 

(2\ ( 3 ) 

The kernels fa 1 a 2 (. x ii X V, fa 1 a 2 a 3 {xi,x 2 ,x 3 ), etc., are posited to have the expansion 

^2) = fol iaa (x 1 ,x 2 ) + e 2 f 2 [2) aia2 (x 1 ,x 2 ) + ••• 
fa%a 3 (x!,x 2 ,x 3 ) = (x x , x 2 , x 3 ) + e 2 ff ^ iaa<l8 (a* , x 2 , x 3 ) + • • • (22) 
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We substitue this into the Schrodinger equation and, by equating coefficients of terms with 
equal number of J's, we get a number of recursion relations. The term corresponding to 
zero powers of J gives a constant term which is a normal ordering term for the Hamiltonian. 



It is taken account of by the normal ordering indicated in the potential energy in (20). The 



coefficient of the term with one power of J vanishes by color contractions. From the term 
with two powers of the current, we find 

2m fgl( Xl ,x 2 ) + 4 [ f£i(x 1 ,x)(nX b (x,y)f£l(y,x 2 ) + V ab 



:r,y 



6 / fifa 2 ab( X ^ X ^^ X ^y)( nU )ab(x,y) +3 / f^' ab (xi , X, y) (Q 1 )^ (x, y, X 2 ) 
x,y Jx,y 



r(3) 



= 
(23) 



For p > 3 the recursion relation is 



mpM~a p + £ n(p + 2 - n)f£U_An% b fl:;r.: 



(n) 



S0\ Ap-n+2) 



n=2 



33- 1 



(p-n+1) 



n=2 
.2 



(p+l)(p + 2) ( p+2 ) ^ i P(P+1) f (p+l) 
afA 1 ' "I o 



(24) 



In these equations, we have used the abbreviations 
(n°) ab (x,y) = 5 ab d y G(x,y) 



(£l l ) abc {x,y,z) 



f*°[6(z-y) + 6(z-x)] G(x,y) 



V ab (x,y) = Sat, d z S(z - x) d z S(z - y) 



(25) 



At the lowest (zeroth) order in e, we have to solve (23) for /q 2 q iO2 0ei, x 2 ); this leads to^j 



m + E a 



(26) 



This agrees with the kernel used in the Gaussian term in (18). Thus, to the lowest order in 



this expansion, we get the same result for the string tension, namely, equation pi). 

In this paper we will outline calculations to the next order, i.e., to order e 2 . For this 
we will need the lowest order results for /( 3 ) and f^. The recursive solution of equations 



L For the holomorphic/antiholomorphic components k,k, we use, k — i(ki + ifa), k — |(fci — ik?). In 



expressions like Ek = \/k 2 + m 2 , k 2 denotes k\ + k 
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(23 25) to order e 2 gives the following lowest order expressions for the cubic and quartic 

ja x a 2 a 3 



kernels. 



24 

ea\a 2 c f b 1 b 2 c 



(2ir) 2 6(k 1 + k 2 + h) g( 3 \k 1 ,k 2 ,k 3 ) 



(27) 



/o axaW* 1 '* 2 ^ 1 '®) = (2tt) 2 J(A: 1 + A; 2 + ?1 + q 2 ) gWfa, k 2 ; q lt q 2 ) (28) 



where 



i {3 Hh,k 2 ,k 3 ) 



64 
16 



k\k 2 (ki - k 2 ) 
E kl + E k2 +E k3 \ {m + E kl ){m + E k2 



+ cycl. perm. 



(29) 



and 



g^\ki,k 2 ]qi,q 2 ) 



E kl + -Efc 2 + E qi + 



</'2 



5 (3) (/ci, fe 2 , - fe) j 1 + 5 (3) (ft- 72- ~<h - 72 ) 



ki + k 2 

(2fci + fc 2 )^i (2fc 2 + fei)^2 



m + E kl 



— g( 3 \ki, k 2 , — fei — fc 2 ) _ 

7i + 72 

These are defined in terms of the Fourier transforms 



m + E k2 
4 



ff (3) (7i,72,-7i -72) 



ki + fc 2 
(2qx + 72) 71 (2g 2 + qi) q 2 



m + E. 



m + E, 



<-ii 



(30) 



fafa 2 a 3 (x 1 ,X 2 ,X 3 ) = J dflfa , . . . , fc 3 ) exp ^ ^ J fafa 2 a 3 ( k l, k 2 , k 3 ) 

/i^U 3 a 4 (^i> ^2, X3,Z4) = y d/i(A?i,...,fc4)exp fiJ^A*^ /iS 2 a 3 a 4 (A;i, k 2 ,k 3l k 4 ), 



where 



d 2 /ci d 2 fc n 
(2^)2 "' (2^)2 



(31) 



(32) 



( 3 ) 

Once again, the value of /q a ^(aJi) x 2 , x 3 ) agrees with previous calculations [2]. Note 
also that .(, ;, (^1, fc 2 ; 7i, 72) as defined in (28 30) is symmertic under independent ex- 
change of the first and second pairs of indices as well as under the simultaneous exchange 
({ai, ki}, {a 2 , k 2 }) <-> ({61, 71}, {b 2 ,q 2 }). We could have certainly made it completely sym- 
metric but we prefer not to do so to keep the notation simple. 

Using the expressions for in (23), the order e 2 -term in /( 2 ) is given by 



/f(7) 



rn 



d 2 k 1 
32tt k 



]( 3) (q,k,-k-q) + 



d 2 k k 
64vr k 



jW(q,k;-q,-k) 



(33) 
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The kernels f^ n \ n > 5, become nontrivial only at the next order. 



Equations(26 33 1 give our recursive solution for the wave function to order e 2 . 



The explicit evaluation of (33 ) presents no difficulties. Two observations greatly simplify 



the task. First, the mass m regulates infrared behavior of the integrals appearing in (33), 
thus enabling the expansion in the external momentum. Second, we are interested only in 
the leading 0(q 2 ) term. This way we find the following analytic result 

/ (2) (?) = + f logf) + = £(1.1308) + ... (34) 

m \ 62 4 2 J 2m 

Comparing this with the coefficient of 0(q 2 ) term in the expansion of the zeroth order kernel 
we conclude that this is equivalent to a —113.08% correction. However, other corrections 
of this order need to be taken into account before we can reach any conclusion. 

3.2 Computation of expectation values 

In carrying out calculations with ^g^o = e 2 , we have to integrate over the currents 

J (and J) with the integration measure for the hermitian WZW theory. This procedure will 
be worked out shortly. It will turn out that the expectation values of products involving only 
the currents J or only the currents J are straightforward to evaluate but those involving 
both J's and J's are more difficult. For this reason, we shall use a slightly different strategy. 

Consider going back to the description in terms of the gauge potentials A, A. The ground 
state wave function is expected to be real, so that the general formula for the expectation 
value of an observable O (which must be gauge-invariant) can be written as 

(O) = JdiM*Z(A,A)0(A,A)MA>A) 

= J dfj,y (A,A)O(A,A)y Q (A,A) (35) 

Now the integrand is a functional of A, A, but we can also write it in terms of J as 



A = M t_ 



+ M^~ l dM^ 

CA 

A = M ] - l dM ] = M t_1 [ ] M f + M ] ~ l dM ] (36) 



$ now becomes a function of and J. Equation (36) shows that we may think of {^4, ^4} 
as the (complex) gauge transform by of {—(tt/ca)J,0}- As argued in [2], we may then 
use the gauge invariance of the wave functions to set to 1, making ^ a function of J. 
(In other words, one may think of our transformation to the J- variables as a complex gauge 
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choice.) Once this is done, equation (36) gives 

(O) = [ dfM(H)e 2cAS —W vfr (J)*o(J) O(J) 



dn(H)e 



2c A S wzm {H) F(J) 



O(J) 



(37) 



The problem is thus reduced to the computation of the correlators of J. (One may wonder 
whether and how it is possible to calculate with J and J. It can be done, in principle; some 
comments about the difficulties of such a procedure are given in the appendix.) 

Now, the current J is constrained in terms of the matrix H by the relation J = 
{cA/ir)dHH ; the fields c/? a in H = exp(t a ip a ) are the unconstrained variables of inte- 
gration. The direct integration of ^q^o ° ver the fields (p a is very involved and tedious. Our 
strategy will be to rewrite the expectation values in terms of a functional integral over an 



(O) 

z 



1 

z 



unconstrained boson field. We start by rewriting (37) as 
O(J) e F ^ 
j dfi(H) e 2cAS —W 

where J a 
1 



dn(H) e 2c * s ^W e -*fC*(8HH-V 



(38) 



c=o 



^2-n/mcA The remaining integral can be evaluated as 

1 



dfi(H)e 



-2c A S wzw (U) 



2c A S wxw (UH)-2c A S wxw (U) 



where C = U 1 dU and we have used the Polyakov-Wiegmann identity 



7T 



TviCdHH- 1 ) = S WZW {UH) - S WZW {U) 



(39) 



(40) 



The expression exp(— 2caS wzw {U)) is the inverse of the chiral Dirac determinant in two 
dimensions. We can therefore represent it in terms of a chiral boson field (in two Euclidean 
dimensions) as 

eM-2c A S wzw (U)) = f [dipdtp] e -S^ B+€ ^ (41) 



The complex boson field (p transforms as the adjoint representation of SU(N). Upon using 

(42) 

S(<p) = J (pdif - F(^2ir/mc A <f~t a <p) 



this back in (38), we find 

(O) = / [dcpdp] e~ s{ ^ 0(^2ir/mc A <ft a V ) 

where 



(43) 
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There is still one correction to be made to this formula. This is because the interactions 
of the chiral boson, represented by F, are such that there is renormalization of the free 



action J ip{d + C)<p. In such a situation, the use of the integral representation (41) for the 



functional integral on the right hand side of (38) is not adequate. At the free level, with no 
interactions coming from e F , we can use (41 ). When interactions are added in, we must ask: 



what representation for the determinant is valid such that we revert to the usual formula 



(41) in the kinematic regime where such interactions vanish? We notice that the F terms 
come with powers of momenta of the currents and vanish as these go to zero. So we must 
use a representation for the determinant which preserves this property. This means that, 



instead of (41), we must use the formula 

exp(-2c A S wzw (U)) = J[d<pd<p] exp - J (p(Z 2 3 + Z x C)y (44) 
The renormalization constants Z\, Z 2 are anyway equal to 1 in the absence of the interactions 



in e , so there is no contradiction with the free formula (41 ), when there are no interactions. 
However, we can also choose Z\, Z 2 to cancel any corrections due to loop integrations arising 
from the interactions e , so that we do recover the free formula when the momenta of the 
currents go to zero. With this more general formula, we get 

S((f) = J {Z 2 CpBip + Zi(pCtp) - F(Z 1 ^/2ir/mc A tpt a <p) (45) 

We could set C to zero at this point, but we have kept it in the above formula to show that 
we can calculate Z\ from the renormalization of the term (pC<p, i.e., vertex renormalization, 
and the same factor should eliminate all divergences from the interactions F. 



Equation (45) is our formula for the action at the chiral boson level. We can now treat 
this as a standard two-dimensional field theory and calculate corrections to the various 
terms. In particular, we are interested in corrections to the term F^ which is quadratic 
in the currents. From the point of view of the chiral boson, this is a four-point vertex. 
The term F^ 2 ' itself does not have any powers of e, so that, to be consistent, we must take 
account of all loops due to this term. This is an additional resummation that has to be 
done. 

As an important example of this resummation, consider the two-point function for the 
currents, which is given, up to constant factors, by (ipt a ip(x) ipt b ip(y)). This may be repre- 
sented diagrammatically as shown in figure 2, where the solid lines represent (<p cp) prop- 
agators and the shaded circle is the vertex corresponding to F^ 2 \ The contribution from 



the free part of the action ( 45 ) , represented by the first term on the right hand side of this 
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V + 



Figure 2: Corrections to the two-point function of currents due to the vertex 



equation, is 



{(pt a <p{x) <pt b <p(y)) 



Tr(t a t 

x ab °A 



j-y 



CA 



tab 



yx 



IT 2 (x — y) 2 



7T 



Jk(x-y) 



k 
k 



Using (26), the contribution of the term with one insertion of the .F^ 2 ) -vertex is 

d 2 k ik(x-y) k 

k 



(i) 



rab c A 



7T 



(2tT 



m) 



m 



The summation of the series of terms shown in figure 2 is thus given by 



(<pt a ip{x) <pt b y{y)) 



?ab c A 



d 2 k 



Jk(x-y) 



k 



rn 



(46) 



(47) 



(48) 



tt J (2vr) 2 k 

The presence of the m/E^ factor improves ultraviolet convergence of integrals and it will 
also suppress the numerical values of various corrections. In any diagram, the vertex cor- 
responding to a current (pt a ip will have such a factor of m/E^. This follows from noticing 
that the current in any diagram has a series of terms (due to insertions of F^) correcting 
it, as shown in figure 3. The summation of these terms gives the result 



ipt a ip(x) 



eff 



<pt a <p(x) - j 


f d 2 k 
1 (2vr) 2 


f d 2 k jk(x- 
J (2^) 2 


-*) 


1 - 


f d k jk( x - 

J (27T) 26 


-*) 


m 
Ek 



Jk(x-z) Ejk_ 



m 



rn 



{(pt a v){z) + 



Ek 



m 



+ 



rn 



Ek 



m 



rn 



+ ... 



(<pt a <p)(z) 



(49) 



For the two-point function for the currents, we must use the corrected current given by this 
equation only at one vertex; otherwise, there will be double-counting. This is similar to the 
case of Schwinger-Dyson equations in, say, electrodynamics, where the vertex corrections 
to the vacuum polarization only apply at one vertex. 
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Figure 3: The effective current vertex 



3.3 The renormalization terms 

We shall now consider briefly the renormalization constants Z\, Z<i- 

The self-energy and vertex corrections for the field (p may be represented by the graphs 
shown in figure 4. Direct computation then shows that the corrected action, to this order, 
may be written as 

4"7tcr f d 2 k I m 



S(ip) = / Z 2 ipdip 



mcA J (2vr) 2 E k + m E k 

+ / Z\(pCtp 



An , f d z k 1 m 

1 - z^r\°R ~ 2 c a, 



mc A ^ 2 ^' J (2n)2E k + mE h \ + "' (50) 
For the moment, we have taken the fields ip to be in an arbitrary representation R. Notice 
that, if we consider only the terms proportional to cr, we have Z\ = Z2, as in electrody- 
namics. However, because of the c^-terms, Z\ and Z2 end up as different constants. Such 
difference between the wave function renormalization and the vertex renormalization is well 
known in nonabelian theories, with exactly the same mismatch as we find here. Reverting 
now to the adjoint representation for the fields, we can identify 

27r f d 2 k 1 m 

Zi « 1 + 



m J (27r) 2 Ek + m E k 
4-7T f d 2 k 1 m 
22 " (2n) 2 E k +mE- k (51) 

The potential divergences are logarithmic; if Z\, Z2 are chosen to cancel them, as above, 
then there are no further divergences. In fact, the corrections to the vertices F^ n \ which is 
what we are interested in, are finite. 



3.4 Evaluating loop corrections 

We have outlined a procedure for carrying out the summation of terms arising from inser- 
tions of the vertex F^ 2 \ which has to be done to all orders. The self-energy and vertex 
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Figure 4: Self energy (upper) and vertex (lower) corrections for the field <p 

corrections of the type shown in figure 4 are compensated by the choice of the renormaliza- 
tion constants. After setting C = 0, there are still factors of Z\ in the vertices F. These are 
compensated by the corresponding insertions of of the i-channel type or vertex type. 
Likewise, there will be self-energy insertions on the propagator lines which are compensated 
by the Z2 factors. The effect of the insertions of F^ is then reduced to terms of the type 
in figures 2 and 3. These can be taken care of by factors of m/Ej~. Thus, we have a fairly 
simple procedure for calculating loop corrections. We write down the loop corrections due 
to the F( 3 \ F^ vertices, etc. Then for the current vertices, we insert factors of m/E^ 
appropriately; each current vertex has such a factor, where k is the difference of the mo- 
menta of the cp and (p legs of the vertex. This counting applies in all cases, except when 
the two-point function of currents arises, in which case there is only one such factor. The 
expressions for the loop corrections, so modified, include the effects of arbitrary number of 
insertions and can be evaluated to give the corrections to various physical quantities. 

4 Corrections to the string tension 

We now turn to the corrections to the string tension which arise from corrections of the 
first type, i.e., from propagator corrections, as shown in figure 1. For this purpose, we must 
calculate, in the (^-language, the corrections to F^ 2 K Actually, for the string tension, the 
low momentum limit of such corrections are adequate. The relevant terms can be easily 
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identified by the Feynman diagrams. To order e 2 , only the F^ 3 ) and vertices are 
relevant, in addition to of course. The basic loop diagrams are as shown in figures 5-9, 
where they are arranged by powers of m/Ek- The action S for the calculations is given by 
equation (45). 

We shall now give the contributions of each diagram to the two-point kernel 
/aiLO^i) x%) = ^aia 2 [/^(<?)]a:i,aj2- Actually, we are interested in the corrections to the string 
tension. For this purpose, it is sufficient to consider contributions to the low-momentum 

(2) 

limit of /q (q), which is given by 

~f (52) 
2m v ' 



/o%) 



m + E g 



Since the corrections generated by the loop diagrams are quite complicated algebraic expres- 
sions, we will give a summary of the low-momentum limits here. Details of the calculations 
will be presented in an appendix. 

First of all, notice that the first correction to from the recursive solution, as obtained 



in (34), gives 

A/ie 2 ]« ^-(1.13082) + ••• (53) 
2m 

The contribution from diagram 1 is given by 

f 

Diagram 1 « —(-0.58118) + ••• (54) 
2m 

This comes from an integral with one power of m/E/.. 

Diagrams 2a and 2b, shown in figure 6, contain two powers of m/Ek- The contribution 
from these diagrams is given as 

Q 2 

Diagram 2a —(-0.47835) + • • • 
2m 

Q 2 

Diagram 2b — (0.20169) + ■■ • (55) 
2m 




Figure 5: Diagram 1 which is the contribution with one factor of m/Ek. 



17 





Figure 6: Diagrams 2a and 2b which are contributions with two powers of m/E^. 

Notice that diagram 2a has the current (pt a ip at one vertex, while from the other two 
vertices we get nonlocal expressions like (p{x)t a (p{y) . Thus we will need to expand this in 



terms of the currents to identify the contribution (55). We have carried out this expansion; 



it is essentially an operator product expansion and will be discussed in appendix B. 

There is one term, namely, diagram 3, which has integrals with (m/E^) 3 . Here too 
we will need an expansion of the product (p{x)t a <p(y) in terms of the currents to pick up 
contributions of the type dJdJ. Once this is done, the contribution of this diagram is seen 
to be 

q 2 

Diagram 3 « —(-0.23569) + • • • (56) 
2m 

There are six diagrams at the level of the (m/E^) 4 , as shown in figures 8 and 9. Of 
these, diagram 4a is zero by the structure of color contractions. The contributions from the 
next four diagrams are as follows. 

Q 2 

Diagram 4b » —(0.02083) + • • • 
2m 

t 

Diagram 4c « —(-0.06893) H 

2 ™ (57) 

Diagram 4d —(-0.01216) H 

2m 

Q 2 

Diagram 4e « —(0.06824) H 

2m 

Diagram 4f has some subtleties and we shall take up its calculation shortly. 




Figure 7: Diagram 3 which is the contribution with three powers of m/E^. 
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Diagram 4c Diagram 4d 



Figure 8: Diagrams with four powers of m/E^ in the integrals. 

At this point, it is useful to take stock of where we are in terms of these corrections. 
It is instructive to look at the result in terms of partial sums of contributions to to a 

(2) 

given order in powers of m/E^. Let A/„ denote the sum of contributions including up to 
n factors of m/Ek m the integrand. We then find Afn = (q 2 /2m)C n , with 

C = 1.13082 
Ci = 0.54964 

(58) 

C 2 = 0.27298 
C 3 = 0.03729 

Notice that the corrections are systematically getting smaller. The first term Co seems 
alarmingly large, but our main point is that it should not be considered in isolation. In 
fact, it is easy to check that it arises from the cubic term in H\ which corresponds to the 
term f abc J c (x)/(x — y) in the operator product expansion of J a (x)J b (y). Diagram 2a, for 
example, has a term which is exactly of this type; if we neglect the factors of m/Ek, it 
would cancel the part of Co arising from exactly. This shows that various contributions 




Diagram 4e Diagram 4f 



Figure 9: Two more diagrams with four powers of m/E^ in the integrals. 



19 



should be considered within certain natural groupings. Our expansion does not put them 
together in such a way from the beginning, but the partial sums are essentially doing this. 
(Individual diagrams give seemingly large values. This is in consonance with [TT]; the form 
of the wave function used there, we suspect, is picking out only the contribution due to 
some of the diagrams we have.) 

The corrected value of the string tension is 



/ _ = e2 /qapb ( 1 + 0.173 + ...) forAjf 



47T 

2 (1 + 0.019 + ---) forA/f 



(59) 



We see that the numerical value of ^/or, up to the diagrams of order (m/Ek) 3 , has ap- 
proximately a 2% correction compared to Q; this correction, however, moves it further 
away from the recent lattice estimates. Of course, we still have to include the diagrams 
with 4 factors oim/E^. However, the main point we want to make at this point is that the 
corrections are small, of the order of a few percent, and that there is a systematic expansion 
scheme which gives a sensible ordering of the various corrections. 

The total contribution from the diagrams 4a to 4e is 

A/£ } _ 4e = ^(0.00798) (60) 
This is a numerically small correction. 

As for diagram 4f, notice that it can be thought of as arising from an effective vertex 
represented by diagram 2a which is then Wick contracted with the vertex. Among the 
diagrams we are considering, this is the only one which involves contraction with an 
vertex. This is why it is special. Now, imagine that we are considering the low momentum 



limit of diagram 2a, namely, A/^ given in equation (55). When further calculations 



are done with the corrected of the form + A/^ , we naturally encounter Wick 
contractions between the fields which go with the two terms. This generates a term which 
is essentially the same as the diagram 4f, but with the loop integration over the new loop 
momenta restricted to low momentum values. (We are first taking the low momentum limit 
and then calculating the subsequent diagram, so the loop integrals in this latter calculation 
have to be cut off at some value. This is why only the low momentum part of the integrations 
occur.) In other words, a part of diagram 4f will be generated in subsequent calculations 
we do in using the effective action for low momentum correlators. Thus, to avoid double- 
counting, only the contribution from the high momentum part of the integration over the 
second loop momentum in diagram 4f should be included as a correction at this stage. 
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Unfortunately, the separation between high and low momentum in the loops is somewhat 
ambiguous. If we calculate diagram 4f with the integrations done only over momenta > m, 
the value is 

A/$ « ^(-0.1037) (61) 
J 2m 



If this value is used, we find 



A/f = £-C A « ^(-.05843) (62) 
2m 2m 



The corrected string tension is then given by 



^ = e y^ (1 - 0.02799 + ••• ) for A/f } (63) 

This is a 2.8% correction with the correct sign, moving it closer to the lattice value. We 
are actually overshooting the lattice value by some amount, but it should be kept in mind 



that the estimate of the last correction in (61) has some ambiguity due to the choice of 



the point at which we should cut off the low momentum integrations. (For the sake of the 
argument, if 2m is used as the cut off, the value of y / o r R differs from Q by —0.00290, about 
one-third of one percent.) Clearly it is important to settle the issue of this cut-off on a 
priori grounds. However, one can see that this issue will not affect the results by more than 
a few percent. An estimate of how bad this ambiguity could be is given by calculating the 

(2) 

integral over all values of the loop momentum. The value of A/^y is then approximately 
(q 2 /2m) x (—0.166). This could decrease the value of y/<jR by another 3%. In other words, 
there is a certain stability to the analysis. The ambiguity of where the low momentum cut 
off should be does not change the value by more than a few percent. 

Can we get the corrected string tension to match the lattice value exactly? First of 
all, for this, it is necessary to analyze the question of the low momentum cut off in the 
integral in diagram 4f more carefully. There are also other corrections such as a number of 
diagrams with 5 and higher number of factors of m/ 'Ej.. There are also corrections of order 
e 4 ; these are expected to remain small because diagrams at this order also have more factors 
of m/Ek- But we are talking about corrections at the one or two percent level and to get 
exact agreement at this level of precision requires consideration of many such contributions. 



5 Summary and Discussion 

The analysis of the Hamiltonian formulation for Yang-Mills theory in 2 + 1 dimensions 
carried out here has helped to place the computation of the wave function and the string 
tension in the context of a systematic expansion. This expansion involves a splitting of 
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the Hamiltonian into Tlo, which is quadratic in the current or functional derivatives with 
respect to it, and Hi, which is cubic in these quantities. The Schrodinger equation is then 
solved recursively in powers of the coupling constant e and the mass parameter m; these two 
quantities are treated as independent at this stage, the relation m = e 2 CA/27r is used only 
at the end. This helps to systematize the expansion parametrically. We have expressed the 
evaluation of expectation values with the vacuum wave function as the functional integral 
of a two-dimensional chiral boson theory. The calculation of low momentum correlators 
can then be visualized as a two step process. First, we define an effective action which 
incorporates the loop corrections in the usual way; this action can be used for calculating 
low momentum correlators as a second step. 

One of the interaction vertices, namely the term which is quadratic in the currents, has 
to be included to all orders in most of the calculations. The effect of this is to introduce 
additional factors of m/Ek (for an appropriate momentum k) in the loop integrals. This 
improves convergence of the integrals and makes their numerical values smaller. As a result, 
it is useful to classify diagrams further by the number of such factors. We have calculated 
the first set of corrections in such a scheme, to first order in the recursive solution for 
(i.e., to first order in e 2 ), including all terms with up to 4 factors of m/E^ in the loop 
corrections which contribute to order e 2 . The result is a correction to the string tension 
which is approximately in the range —0.3% to —2.8% in ^/or, compared to the value given 
in This makes the result consistent with the recent lattice estimates. 

Regarding theoretical uncertainty, we note that there is an ambiguity in one of the 
integrals, as to precisely where the separation between what is identified as "low momentum" 
and what is identified as "high momentum" is to be imposed. This needs to be clarified, 
even though we see that the variance due to this can be no more than a percent or two. 
There are also corrections involving 5 or more powers oim/E^. There are also higher terms 
(of order e 4 ) in the recursion and corresponding loop corrections. All these corrections are 
expected to be small, based on the fact that there are more powers of m/Ek, but when 
talking about corrections of the order of a percent, they can have a significant effect. 

The true value of our analysis is to show that corrections can be systematized and 
remain under control and are numerically small. In this sense, we have demonstrated the 
feasibility of systematic calculations in our Hamiltonian approach to Yang-Mills theory. The 
final value of the string tension, to the order we have calculated, is also consistent with the 
recent lattice estimates. 

This work was supported in part by the National Science Foundation grants PHY- 
0555620 and PHY-0758008 and by PSC-CUNY grants. 
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APPENDIX A 



Our calculations have all been done in terms of correlators of J. Once the wave function has 
been calculated in terms of J, Vl/* naturally involves J and one may ask whether it is possible 
to do calculations using 'I'* 'I'. In principle, this is definitely possible. Since the complete 
wave function ^o(J) is expected to be holomorphically invariant, we can re-express ^q(J) 
back in terms of J, and we only have to evaluate correlators of J. The difficulty, however, 
is that the separation of the Hamiltonian into T~Cq and H\ is not holomorphically invariant 
and since we are only calculating to a certain order in e 2 , we do have to consider the 
issue of expressing J in terms of J. (The fact that this separation is not consistent with 
holomorphic invariance is not a problem in itself; it preserves the invariance in an order- 
by-order fashion. Recall that a similar separation is standard in gauge theories whereby 
one loses the full gauge invariance at a particular order. However, gauge invariance is 
preserved, in a consistent fashion; i.e., up to terms of higher order. The Ward identities, 
likewise, connect vertices and Green's functions of different orders. A similar situation is 
expected in our case.) 

Now, from the definition of the currents, dJ a = dJ b H ba , where H ba = 2Ti(t b Ht a H- 1 ) 
is the field H as an adjoint matrix. One way to proceed would be to rewrite J a in terms of 
J a as J a = ±(dJ b H ba ) and further express H in terms of J as an infinite series, ^q^o then 
becomes a function of J's and expectation values can be evaluated. The difficulty with this 
approach is that one may need to take account of the full infinite series to get meaningful 
results. As an example, recall the calculation of the two-point function of the group- valued 
field g in a unitary WZW model. We can write g = exp(it a 6 a ) and expand in powers 
of a and calculate the correlators. To the lowest order this gives a logarithmic function, 
(g( x )g(y)} ~ (8 a (%)& b (y)} ~ 5 ab log(x — y). However, the whole series can be summed up to 
give a power law result which incorporates the anomalous dimension of g correctly. (The 
summation is usually done by more efficient means such as the Polyakov-Wiegmann identity 
or the Knizhnik-Zamolodchikov equation, but that is a separate issue.) In the case of our 
field H, we expect a similar situation. We have checked that at least some of the lower 
order terms lead to logarithms which are likely to be spurious. The full resummation is, 
however, beyond what can be done at present. Therefore, the method used in subsection 
3.2 is significantly better. 

If calculations are done, taking into account the relationship between J and J in full, 
then we expect that the two ways of calculating, namely, using ^q^q with J transformed 
back to J's as above, or writing the expectation value in terms of J by setting to 1 as 
done in subsection 3.2, will yield the same result. 
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APPENDIX B 



In this appendix we present details of the computation of the diagrams shown in Figures 
51(8} We start by collecting various expressions relevant for this discussion. 



We consider the two-dimensional chiral boson theory ( 43 ) with interaction term 



2 

F((pT a ip) = F^^T a ^) + ^F ( - 3 \^T a ip) + jF^(^T a ip)+ ... (Bl) 
where interaction vertices F^ 2 \ F^ 3 \ are given by 

F (2) = Jlf Mq) ( ^ TV) } ^T^U (B2) 

mcA J 

F (3) = ( — J / dv(h,k 2 ,k3) (ipT a ^)- kl {CpT a *<p)„ k2 ((pT a °v)-k 3 



ttica 



Xf0 3) ai a 2 a 3 (kl,k 2 ,k 3 ) (B3) 

F(4)= i^t) I Mh,k 2 , qi , q2 ) (^r a v)- fel (^ a ^)_ fe (^T 6 ^)_ 3l (^r 6 ^)_ 92 (B4) 



with kernels / (2) , / (3) , / (4) as in eqs. (26 30), and 



M<l) = -^ 2 : d ^ kuk2) = ^2^^' etC - (B5) 
For the Wick contraction of the chiral boson fields and currents we have 

(^ a (r) tp b (s)) = S ab (2ir) 2 5(r + s) — 

ir 

({<pT a tp) k (<pl*<p) p ) = 5 ab {2n)H{k + p) (B6) 

Diagram 1 

Diagram 1 in Figure [5] is quite easy to evaluate. Contraction of one pair of bosonic currents 
in F^\(pT(p) leads to the following expression 

Diagram 1 = ^- J d^q) (0TV)-, {CpT a ip) q g(%, k; -q, -k)\ (B7) 



The only less trivial step in deriving (|B7|) is that combinatorial factor is = 4 instead of 

2 

the quartic kernel ffl^^ki, k 2 ; q\,q 2 ) 



the naive expectation = 6. This is so because we chose not to symmetrize completely 
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Comparing (B7) with F^ 2 \(pTip) from eq.(B2) we immediately conclude that correction 



to quadratic kernel f^ 2 \q) due to this diagram is 



A/«(g) = J 



d 2 k k m 
64^ k Ei 



9 



W(q,k; -q-k) 



(B8) 



The leading, 0(q 2 ), term in the low-momentum expansion of (B8) can be computed ana- 
lytically; we find 



2m 



(-0.581178) + ... (B9) 



Diagram 2b 



Of the two diagrams with two factors of m/E k shown in Figure [6] , diagram 2b is straight- 
forward, so we will evaluate it first. The Wick contraction of two pairs of bosonic currents 
gives 



Diagram 2b 



2tt 
mcA 



d^q) ($>TV)_, 



d 2 k k(k + q) m 



128tt k(k + q) E k E g+k 



P(q,k,-q-k)] (BIO) 



Again, comparing this expression with (B2) we conclude that correction to f^ 2 \q) due to 
Diagram 2b is 



A/< 2 >(g) 



d 2 k k(k + q) m 



128tt k(k + q) E k E q+k V 



p\ q ,k,-q-k) 



(Bll) 



The leading, 0(q ), term in the expansion of (Bll) can be calculated analytically as 



= I; (-f + £ tag*-'- ? log* i + .. 



2m 



(0.201688) + 



(B12) 



Diagram 2a 



Diagram 2a provides the first nontrivial example of diagrams we have to deal with. First 
of all, notice that it is based on a single cubic vertex and, as such, this diagram does 
not have a direct analogue in standard perturbative QCD. Second, while it is elementary 
to write an explicit expression which corresponds to this diagram, i.e., 

27T f 

Diagram 2a = / dfi(q, r) (^T a v7)_„te(r)T a 93(-r+g)+^(-r+g)r a ^(r)]n(r, q) (B13) 

mc A J 



with 



U(r,q) 



d 2 k g( 3 \k,q,—k — q) m 2 
64^ ~ 



r + k E k E k+q ' 



(B14) 
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the interpretation of this expression requires some care. The problem is that the term which 
is bilinear in (p and <p, namely, 

d 2 r 



1 = 



(27T) 



[¥{r)T a <p{-r + q)+ <p(-r + q)T a <p(r)\ U(r, < 



(B15) 



is not simply proportional to a current {(pT a ip) q and therefore (B13) cannot be directly 
compared to (B2). This is because H(r,q) depends on r in a nontrivial way We note 
that n(r, q) admits a power series expansion in r, which means that, in addition to the 



current ({pT a ip) q , the expansion of (B15) also contains other local operators of the type 
{d n d m (pT a tp) q and {(pT a d n d m (p) q . But, in our two-dimensional theory, any such local op- 
erator can (in principle) be written entirely in terms of currents (and products of currents) 



only. In other words, Diagram 2a and equation (B13) generate not only corrections to the 



quadratic kernel f^ 2 \q) but also to the cubic and higher kernels as well. In this paper 
we are only interested in the leading 0(q 2 ) correction to f i - 2 \q) which corresponds to the 



q (ipT a (p) q -teTm in the operator product expansion of (B15). Therefore we will concentrate 



on extracting this term from (B15). 

The integral Il(r, q) admits a power series expansion in both arguments, so we may write 



II(r, q) = Aq 2 + Bqr + Cr 2 + ... 



(B16) 



where A, B and C are some constant coefficients^] and the ellipsis stands for higher order 
terms (like q 2 fq, f 3 r etc.) in the expansion. Note that these higher order terms have at least 
one power of holomorphic momentum q or r. Therefore, simple power counting suggests 



that such terms will not contribute to the q 2 J^-term in the OPE of (B15). 



We will consider the contribution (to (B15 )) of each of the three terms on the right hand 



side of (B16) separately. The A-term is trivially evaluated as 

<hr [^(r)TV(-r + q) + (p(-r + q)T a cp{r)} Aq 2 = 2 A q 2 J a q 



(B17) 



(2tt)2 

Similarly, the B-term can be evaluated as 

^ W)T a v{-r + q) + (p{-r + g)TV(r)] Bqr = Bq [(«^T<V) 9 + {CpT a d<p) q ] 



B q 2 J a Q 



(B18) 



2 Imagine for a moment that we may forget about the r-dependence of IT(r, q) and take, for example, 
II(0,q) instead. In this case (B15I nicely factorizes into 2 11(0, q) (<fT a ip) q , and we see that ( B15 1 would be 



proportional to a current (ipT a ip) q 



Numerical values of these constants can certainly be found from (B14l. However as will be seen shortly, 



we do not really need to do that. 
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The term involving Cr 2 requires more care. Notice that by a change of variables r = —k+^q 
in the first operator product and f = k + in the second one we may write 



/ 



d 2 r 

[ip(r)T a v{-r + q) + (p(-r + q)T a <p{r)} f 2 



(27T) 2 
2 



= j w? ^ { ~ k + q/2)TMh + q/2) + q/2)2 +Ch ~ 9_/2)2] (B19) 

For the composite operator (p{—k + q/2)T a ip{k + q/2), q is the total momentum and 2k is 
the relative momentum of the two fields (p and <p. Notice that powers of k generate terms 
which are asymmetric between ip and (p in coordinate space, such as (p d T a ip; such operator 
structures are not of interest to us at present. The term of interest corresponds to k = 0, 
of the form lq 2 ip(q/2)T a ip(q/2). In other words, we may approximate 



2 

/2 



C [(B 2 (pT«y) q + {(pT a B\) q ] = %fJ\ (B20) 



2 

Since 

2 Aq 2 + B q 2 + \Cq 2 = 2{Aq 2 + Bqr + Cf 2 ) 



_ 1 _ 
r= 5 9 



(B21) 



we may simplify the contribution (B13) of Diagram 2a as 



2ir f 

Diagram 2a = / dfi(q)(ipT a ip)JipT a ip) q 2 U(q/2, q) + • • • (B22) 

mc A J 



which in turn may be directly compared to ( B2 ) to give 



J 327T g/2 + k E k E k+q 

Computation of the leading 0(q 2 ) term in the above integral is straightforward and can be 
done analytically to obtain 

A/ (2) (^) = t~ (~\ + 9 log 2 - 6 log 3^ + . . . = f - (-0.478349) + . . . (B24) 
2m \ o J 2m 



Diagram 3 

Computation of Diagram 3 is quite similar to Diagram 2a which we have discussed at some 
length in the previous section. Straightforward contractions dictated by Feynman graph in 
Figure 7 lead to the following expression 

Diagram 3 = — / (#TV)_, [v(r)T a ^(-r+q)+ip(^+q)T a ip(r)]U(r, q) (B25) 

mc A J (2tt) 4 
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where now by n(r, q) we mean the following 



d 2 kd 2 l g( 4 \q,k-l,-q-k-l) 



m 



3 



""'•' /) ' 2(16vr)2 ( f + £)(r + fc + I) E^E^^ (B26) 



We see that equation (B25) is essentially identical to (B13) and therefore all arguments 
from previous section apply to Diagram 3 as well. In particular, the correction to quadratic 
kernel f^ 2 \q) is 

Af (2^ o nY /o A f d 2 kd 2 l gW(q,k;l,-q-k-l) m 3 

Af ( >(q) = 2Il(q 2,q) = / 2 , r , n F F r B27 

7 (167r) 2 (g/2 + k)(q/2 + k + I) E k EiE k+i+q 

The only difference from the previously considered example is that now we have to evaluate 

a two-loop integral which is difficult to do analytically. However, it is quite easy to evaluate 



(B27) numerically and we find 

a 2 



Af (Z, (q) = —(-0.23569 ± 0.00001) + . . . (B28) 
2m 



Diagram 4a 



When evaluating Diagram 4a we will find that it is proportional to the following color factor 

ja!a 2 a 3 jb 1 b 2 b 3 r^^-pa 2 rpb 2 rpa 3 rpb 3 ^ _ 1 jaia 2 a 3 jbib 2 b 3 r^t^rpa 2 rpb 2 ^rpa 3 rpb 3 j 

_|_p^ a 2 ipb 2 ^^rpa 3 T b 3 y\ 

= (B29) 

since the symmetric tensor d abc = Tr({T a , T b }T c ) = for the adjoint representation. There- 
fore we conclude that, by virtue of the color contractions, 

Diagram 4a = (B30) 

Diagram 4b 

Once we do all bosonic line contractions as prescribed by Feynman graph in Figure 8, we 
obtain the following expression 

2vr f d 2 q 



mcA J (27r) 2 



* 3 ^ <wt ,,, p; » m h.*.-»-*)» m h.'.-»-o 



64m J r(r — q)(r + k)(r + I) 

(B31) 

Notice that the expression in square brackets is a three-loop integral that we have to eval- 
uate. One of the loop integrals can be done analytically, namely, 



d 2 r 1 1 



(27r) 2 r(r — q)(r + k){r + I) ir 



k I 
+ 



k{k + q){l-k) l{l + q){k-l) q{k + q){l + q 

(B32) 
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Using this result as well as inserting 4 appropriate m/E factors into (B31) we end up with 

2vr f d 2 q 



Diagram 4a 



d 2 kd 2 i 



m 



mcA J (27r) 2 

g {3 \q,k,-q-k) g W( q ,l,-q-l) 



(32vr) 2 E k EiE k+q E l+q 
k 



k{k + q){l- k) 



+ 



l(l + q)(k-l) q{k + q){l + q) 

(B33) 

And this translates straightforwardly into the following correction to quadratic kernel f^ 2 \ 



AfW(q) 



d 2 kd 2 l 



(32n 



- gW(q,k,-q-k) g W(q,l,-q-l) 



m 



k 



+ 



I 



EkEiEk+ q Ei +q 

q 



(B34) 



k{k + q){l-k) l(l + q)(k-l) q{k + q){l + q) 
This (two-loop) integral can be easily evaluated numerically and we find 



Af^tq) = — (0.020828 ± 0.000002) + 
2m 



(B35) 



Diagram 4c 



For Diagram^] 4c we easily obtain 
2tt 



Diagram 4c 



mcA 



dfi(q, r, s) [(p{r)T a (p(-r - q) + tp(-r - q)T a <f(r)} U(r, s, q) 

x [<p(s)T a <p(-s + q) + <p(-8 + q)T a ip(s)) (B36) 



where n(r, s, q) stands for 



U(r,s,q) 



d 2 kd 2 l gW(k,q-k;l,-q-l) 



m 



Ek Ei E k -q Ei +q 



(B37) 



(64vr) 2 (k + r){l + s) 

Straightforward generalization of the analysis presented for Diagram 2a suggests that we 
should evaluate U(r,s,q) at r = — q/2 and s = q/2. In other words we may write the 



relevant term of ( B36 ) as 

2vr f d 2 



Diagram 4c = — / {yT a V )„ q {(pT a y) q 4 U(-q/2, q/2, q) + ■ ■ ■ (B38) 

niCA J {^) 



4 We consider Diagrams 4c and 4d as being different because we chose not to symmetrize completely the 
quartic vertex . Had we chosen to work with a completely symmeterized form of F^' , both diagrams 
would have been treated identically. 
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from which we immediately conclude that 



d 2 kd 2 l g^\k,q-k;l,-q-l) 



rn 



l 



Th integral in ( B39 ) can be evaluated numerically and we find 



AfW(q) = — (-0.06893 ±0.00002) + ... (B40) 
2m 



Diagram 4d 



The mathematical expression for diagram 4d splits into two parts involving symmetrized 
and antisymmetrized products of the color matrices. Thus 

Diagram 4d = [Diagram 4d]s + [Diagram 4d]^ (B41) 

where 

[Diagram 4d] s = ~ (^^) / d ^ k ^ k<2 > 51 ' q2 > r ' s "> -^iLW^ 1 ' k ^ q2 "> 
x [^{r){T a \T^} ( — L- + — ^ cp(-r - A* - ft)] 

x &( s ){T a \T b >} ( — L- + — ^ p(- a - - <&)] (B42) 
\s + k 2 s + q 2 J 

The color structure here is such that we will not generate a term of the current-current 
form. If we expand the symmetrized product of the adjoint matrices, say, {T ai ,T bl } in 
terms of a basis of hermitian matrices, the term proportional to T c will have the coefficient 
d aiblC = Tr{T ai , T bl }T c , which is zero for the adjoint representation. Thus the contribution 
involving the current-current form comes from the second part of the diagram with the 
antisymmetrized products. It is given by 

C 2 / 27T \ 2 f 

[Diagram 4d] A = - — ( J / dfi(ki, k 2 , qi, 92, r, s) /^.^(fci, fa; qi, 92) 

x [0(r)[T»STfc] ( — L- - -^r) <p(-r - h - qi )} 
\r + ki r + qij 

x [<p(s)[T a >,T b *] ( zt~t~ - j^zr) <p(rs -fa- q 2 )] (B43) 
\s + k 2 s + q 2 J 

This expression can be further simplified as 

27T f 

[Diagram 4^ = / du(q,r,s) [(p(r)T a ip(-r-q)+(p(-r-q)T a ip(r)]Il(r,s,q) 

mc A J 

x [^)T>(-s + q) + <p{-8 + q)T a tp{s)\ (B44) 
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Figure 10: "Effective" quartic vertex, 
where, we have once again, 

d 2 kd 2 l g^(k,l;q-k;,-q-l) 



in 



A 



n fc«>9) = T^TvI /7" — \ tT i' -\ T-i t-i T-i (B45) 



(64^) 2 (A; + f)(Z + s) EkEtEk^Ei 



1 



The rest of the analysis is similar to the case of diagram 4c, giving the contribution 



d 2 kd 2 l g( 4 \k,l;q-k,-q-l) 



in 



A 



AfW(q) = 4U(-q/2,q/2,q) = / 2 ,f ''l* j-'' = = (B46) 

(32vr) 2 (k - q/2)(l + q/2) E k EiE k _ q E i+q 



Numerical evaluation of the integral gives 

^2 



A/ (2) (g) = — (-0.01216 ±0.00007) + ... (B47) 
2m 



Diagram 4e 



The simplest way to compute Diagram 4e is to introduce an "effective" quartic vertex 



g^/j (fa,k 2 ;qi,q 2 ) = g w (fa,k 2 , -fa - k 2 ) 



fa + k 2 1 



fa + k 2 E kl . 



9 {3) (qi,q2,-qi~q2) (B48) 



as shown on Figure 10 and to notice that Diagram 4e is quite similar to Diagram 3. 



Therefore, we may simply re-use previously derived equations (B25 B27) by replacing g^ 
with g^ff In this way we immediately find 

2ir f 

Diagram 4e = / dfi(q, r) ({pT a ip)„„ [fi(r)T a ip(-r + q) + tp(-r + q)T a ip(r)} Il(r, q) 

mc A J ' 

(B49) 
(B50) 



with 



U(r,q) 



d 2 kd 2 l 9, 



(4) , 

eff [ 



,k;l,—q — k — I) 



m 



2(16tt) 2 (f + k)(f + k + l) EkEiEk+i+g 
As usual, correction to f^ 2 \q) is 



AfW(q)= 2I%/2,, 



m 



d 2 kd 2 l g^ f (q,k;l,-q-k-l) 

(16vr) 2 (q/2 + k)(q/2 + k + l) E k Ei E k+l+q 



(B51) 
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Numerical evaluation of this integral presents no difficulties and we find 

ri 2 



Af {Z) (q) = ^(-0.06824 ± 0.00002) + . . . (B52) 

_■ III 



Diagram 4/ 



Diagram 4f has the structure of the product of an effective vertex corresponding to diagram 
2a and the vertex with two Wick contractions connecting them. It can also be directly 
written down from the interaction terms. We get 

2ir f 

Diagram 4f = / dfi(q, r) (<pT a (p)_ q [<f{r)T a ip(-r + q) + <p(-r + q)T a p{r)} U(r, q) 

mcA J 

(B53) 

where 

f Wi I g$)(k,q,-k-q) ( 1_ 



256vr 2 (f+l)(f+l-q) \f+l + k q + k-f-l 

x- ^ (B54) 

(m + E t )EiE k E k+q 

We can now follow a similar line of reasoning as we did for diagrams 2a, 3, 4c-4e and 
conclude that 



Af^ 2 \q) = 2U(q/2,q) 



I 



d 2 kd 2 l (3)n P , 

gW(k, q, -k - q) -=- — — + 



128vr 2 * ^ (P-(g/2) 2 ) \{q/2)+l + k (q/2) + k-l, 

(B55) 



m 2 



(m + E^EtEkEk+g 
The numerical evaluation of this integral gives 

a 2 

Af( 2 \q) = ^-(-0.1666 ± 0.0002) (B56) 

This diagram, as explained in text, has some subtleties. We have argued that, to avoid 
double counting of a part of this diagram, we have to restrict the integration over the second 
loop momentum k to values above the low momentum cut off. Taking the value of the cut 
off as m and 2m in turn, we find 

, q 2 { (-0.1037 ±0.0003) for cutoff = m 

Af (2 \q) = -f- { K ' B57 

2m 1 (-0.051104 ± 0.00019) for cutoff = 2m 
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